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The Carnot statement of the second law of thermodynamics poses an upper limit on the effi¬ 
ciency of all heat engines. Recently, it has been studied whether generic quantum features such as 
coherence and quantum entanglement could allow for quantum devices with efficiencies larger than 
the Carnot efficiency. The present study shows that this is not permitted by the laws of thermo¬ 
dynamics - independent of the model. We will show that rather the definition of heat has to be 
modified to account for the thermodynamic cost for maintaining non-Gibbsian equilibrium states. 
Our theoretical findings are illustrated for two experimentally relevant examples. 
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I. INTRODUCTION 

Harnessing energy stored in inaccessible forms such as 
heat or chemical energy and transforming it into useful 
work is one of the most important, technological achieve¬ 
ments. Nevertheless, the underlying principles and ulti¬ 
mate limitations imposed by quantum mechanics on such 
thermodynamic processes are still an active field of re¬ 
search [1, 2]. Modern studies range from implementa¬ 
tions of quantum heat engines in ion traps [3] over ther¬ 
modynamic cycles in optomechanical systems [4, 5] to the 
description of the principles of photosynthesis as photo- 
Carnot engines [6, 7]. The natural question arises how 
generic quantum features such as coherence and entangle¬ 
ment affect classical formulations of the thermodynamic 
axioms. This problem has been studied from many differ¬ 
ent perspectives [8] including stochastic thermodynam¬ 
ics [9-11] and information theory [12-14]. However, a 
conclusive consensus appears still to be lacking. 

In particular, it has been studied whether quantum 
correlations could be harnessed, and whether quantum 
devices could operate with efficiencies larger than the 
Carnot efficiency [15-18] - therefore demanding a refor¬ 
mulation of the Carnot statement of the second law of 
thermodynamics. 

The Carnot statements of the second law of thermo¬ 
dynamics declares [19] 

No engine operating between two heat reservoirs can be 
more efficient than a Carnot engine operating between 
those same reservoirs. 

Recent studies, however, raised the question whether 
quantum effects such as coherence and entanglement 
could provide means to break the limit posed by the 
Carnot efficiency [14-16, 20]. To this end, a variety 
of theoretical and experimental setups have been devel¬ 
oped [3, 4, 21]. 
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FIG. 1. Illustration of a generic quantum heat engine: A 
quantum particle inside in a quantum piston consisting, e.g., 
of an optical cavity. 


Thermal equilibrium states of classical systems are 
universally described by the Boltzmann-Gibbs distribu¬ 
tion [22, 23]. It can be easily shown that the Carnot 
statement is only a consequence of this universality of 
equilibrium states [22]. The situation is dramatically 
different for thermally open quantum systems. Gener- 
ically, quantum systems, which are not only ultra-weakly 
coupled to their environment, do not relax into Gibbs 
states [24, 25]. This can be seen most clearly for the 
analytically solvable case of quantum Brownian motion 
[26]. It has been shown that the quantum correlations 
between system and environment prevent relaxation into 
the Gibbs equilibrium states [27]. Therefore, such non- 
Gibbsian equilibrium states are not fully “thermalized” 
and contain additional information, encoded in quantum 
coherence and entanglement [28, 29]. 

It is only natural to ask whether a quantum heat en- 
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gine such as in Fig. 1 could be realized that utilizes this 
“extra” information and thus constitutes a device oper¬ 
ating with an efficiency larger than predicted by Carnot 
[15]. The present study will elucidate that this is not 
permitted by the laws of thermodynamics. In particu¬ 
lar, we will argue that there is a specific thermodynamic 
price that has to be to paid to maintain quantum corre¬ 
lations and thus to prevent the system from relaxing into 
a Gibbs state. How to properly modify the definition 
of heat has been recently studied with great intensity for 
various model systems [5, 30-33]. However, to the best of 
our knowledge, the only rigorous and model independent 
distinction has been previously discussed in the context 
of the heat capacity [34] of open quantum systems (see 
Eq. (23) of Ref. [34]). 

In the present work we further develop this notion of 
quantum heat. As a main result we will see that the 
classical Carnot statement remains unchanged for open 
quantum systems with arbitrary coupling to their envi¬ 
ronment. More specifically, we will show 

No quantum heat engine operating in quasistatie Carnot 
cyeles ean harness quantum eorrelations. 

This insight has far-reaching consequences for all areas 
of engineering at the nanoscale. Nanoengines perform¬ 
ing beyond the Carnot limit necessarily operate far from 
thermal equilibrium [35]. To the best of our knowledge, 
however, this has only been proven rigorously, if at least 
the initial state is Gibbsian [36-38] . 


II. QUASISTATIC PROCESSES 

Before we analyze the quantum Carnot cycle, let us es¬ 
tablish an important concept, first. Consider a quantum 
system with Hamiltonian H(ujt), where Wt is an exter¬ 
nal control parameter. Then the dynamics of the system 
is governed by the Liouville type equation p = 
where the superoperator reflects both the unitary dy¬ 
namics generated by H and the non-unitary contribution 
induced by the interaction with the environment. We 
further have to assume that the equation for the steady 
state, Cujtip) = 0, has a unique solution [39]. The clas¬ 
sical Carnot statement is formulated for cycles of qua¬ 
sistatic (infinitely slow) processes, i.e., successions of sta¬ 
tionary states [40]. In complete analogy to the classical 
theory quasistatic processes are the only processes con¬ 
sidered in the present analysis. 


III. THERMODYNAMICS OF GIBBS 
EQUILIBRIUM STATES 

We begin with a brief review of fundamental thermo¬ 
dynamic concepts for Gibbs equilibrium states, 

p = exp (—/3H)/Z, where Z = tr {exp {—/3H)} . (1) 


Here /3 is the inverse temperature of the environment, 
[3 = 1/T, and we work in units for which the Boltzmann 
constant is unity. The thermodynamic entropy is then 
given by the Gibbs entropy [22], S = — trjplogp} = 
(3 {E — F), where E = tc {pH{ is the internal energy of 
the system, and F = —T log Z denotes the Helmholtz 
free energy. For isothermal, quasistatic processes, /3 = 0, 
the change of thermodynamic entropy dS” becomes 

dS = (3 {tv{5pF[} + {tv{p5H} — dF)) = ptr {SpH} , 

( 2 ) 

where (5 denotes an infinitesimal change. Therefore, two 
forms of energy can be identified [1]: heat is the change 
of internal energy associated with a change of entropy; 
work is the change of internal energy due to the change of 
an extensive parameter, i.e., change of the Hamiltonian 
of the system. We have, 

dE = SQ + SW = tr [SpH} + tT {pSH} . (3) 

The first law of thermodynamics (3) is a universally valid 
expression of the conservation of energy. However, the 
identification of heat SQ, and work SW (3) is consistent 
with the second law of thermodynamics for quasistatic 
processes (2) if, and as will shortly see, only if p is a. 
Gibbs state (1). 

It is worth emphasizing that for isothermal, quasistatic 
processes we have, 

dS = f3SQ and dF = SW , (4) 

for which the first law of thermodynamics takes the form 

dE = TdS + dF. (5) 

In this particular formulation it becomes apparent that 
changes of the internal energy dE can be separated into 
“useful” work dF and an additional contribution, T dS”, 
reflecting the entropic eost of the process. 

IV. THERMODYNAMICS OF NON-GIBBSIAN 
EQUILIBRIUM STATES 

For systems, whose equilibrium states are not de¬ 
scribed by a Boltzmann-Gibbs distribution (1), the iden¬ 
tification of heat only with changes of the state of the 
system (3) is no longer possible [34]. Mathematically 
similar situations have been studied for classical systems 
under non-conservative forcing [41-43]. Such systems re¬ 
lax into so-called nonequilibrium stationary states, and 
it has been recognized that not all heat absorbed by the 
system accounts for the entropic cost [42]. Some con¬ 
tribution to the total heat, coined housekeeping energy 
(or heat) SQuk [41-43], fulfills the sole purpose of pre¬ 
venting the system from relaxing into the thermal Gibbs 
state. The concept of quantum housekeeping heat has 
been analyzed carefully in Refs. [44, 45]. 

For generic quantum systems, whose thermal equilib¬ 
rium states are non-Gibbsian, the situation is mathemat¬ 
ically analogous. However, we stress that in the present 
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context we are interested in quantum systems in non- 
Gibbsian equilibrium states, and not in generic nonequi¬ 
librium situations. The deviation from the Gibbsian 
equilibrium is only due to the interaction and correla¬ 
tions between system and environment. In particular, 
there is no continuous supply of energy from the envi¬ 
ronment into the system. For instance, it has been seen 
explicitly in the context of quantum Brownian motion 
[27] that system and environment are generically entan¬ 
gled. In such situations the reduced equilibrium state, u, 
of the system only can be written as [24], 

cr = exp (-/3 (i7 -I- A))/Za (6) 

where H is the reduced Hamiltonian and the modified 
partition function. A similar situation is encountered in 
classical systems with non-negligible interacting energy 
between a system of interest and its environment. For 
such classical systems it has been shown that H + A can 
be interpreted as a potential of mean force [46], and that 
the identification of thermodynamic work is subtle [46]. 

Note, however, that generically the physical situation 
is even more involved in the quantum case. Whereas 
A for classical systems only includes contributions from 
interaction energies, solvation energies and classical cor¬ 
relations, for quantum systems A is also governed by 
quantum correlations. This means that even for situa¬ 
tions for which the surface terms such as the interaction 
energy are vanishingly small, the purely quantum part of 
A can not necessarily be neglected [24]. 

Therefore, to formulate thermodynamics consistently 
the energetic back action due to the correlation of sys¬ 
tem and environment has always to be considered care¬ 
fully [27, 34]. During quasistatic processes parts of the 
energy exchanged with the environment are not related 
to a change of the thermodynamic entropy of the system, 
but rather constitute the energetic price to maintain the 
non-Gibbsian state, i.e., coherence and correlations be¬ 
tween system and environment. 

In complete analogy to stochastic thermodynamics we 
identify the thermodynamic entropy with the von Neu¬ 
mann entropy [42, 44, 47-49]. With a being the station¬ 
ary state, we can write 


where we identified the total heat as SQtot = ti{SaH} 
and the correlation part as SQc = dJ" — tr {aSH}. 

The total heat exchanged with the environment has 
two contributions. The correlation heat is the energetic 
price that has to be paid to maintain coherence and quan¬ 
tum correlations. The excess heat (5(5ex is the only con¬ 
tribution that is associated with the entropic cost, 

d'H = l3SQe^, and = SQtot - SQc ■ (9) 

Notice that SQex is mathematically analogous to the ex¬ 
cess heat for classical systems under non-conservative 
driving [42]. 

Accordingly, the first law of thermodynamics takes the 
form 

d£ = JWex + <5Qex (10) 

where SWe^ = 5W + SQc is the excess work [47], which 
reduces in the classical limit to the notion analyzed in 
Ref. [46]. Finally, Eq. (5) generalizes for isothermal, qua¬ 
sistatic processes in generic quantum systems to 

d£ = Tdn + dT. (11) 

In the remainder of the present analysis we will show 
how the universal Garnot statement follows from these 
generalized thermodynamic relations. 

V. UNIVERSAL EFFICIENCY OF QUANTUM 
CARNOT ENGINES 

Imagine a generic quantum system that operates be¬ 
tween two heat reservoirs with hot, Thotj and cold, Tcoidj 
temperatures, respectively. Then, the Carnot cycle con¬ 
sists of two isothermal processes during which the sys¬ 
tems absorbs/exhausts heat and two thermodynamically 
adiabatic, i.e., isentropic strokes during which the exten¬ 
sive control parameter uj is varied. 

During the first isothermal stroke, the system is put 
into contact with the hot reservoir. As a result, the excess 
heat Qex,i is absorbed at temperature Thot and excess 
work IFex,i is performed. 


T-L = —tr {(T log cr} -I- (tr {tr log p} — tr {cr log p}) 
= P{£-{F + TD{a\\p)))=P{£-F), 


( 7 ) 


kFex.i = •7 ^(w 2, Thot) — •7-'(a;i, Thot) 

Qex.l = Thot ('H(w2,Thot) — 'H(wi, Thot))- 


where 8 = tr{aH} is the internal energy of the system, 
and SF = F+TD{a\\p) is the information free energy [47]. 
Here, D (ct 11 p) = tr {cr (log a — log p)} is the quantum rel¬ 
ative entropy [50]. Note that it has been shown that F 
is the only thermodynamically consistent definition of a 
free energy for non-Gibbsian states [47, 51]. 

As before (2) we now consider isothermal, quasistatic 
processes, for which the infinitesimal change of the en¬ 
tropy reads 

diH = fS [tv {SaFl} + [tv {fjSFl} — dF)] 

= H5Qtot-SQc). 


Next, during the isentropic stroke, the system performs 
work Wex ,2 and no excess heat is exchanged with the 
reservoir, AT-L — 0. Therefore, the temperature of the 
engine drops from Thot to Tcoidj 

Wex ,2 = A 8 = 8 (uj 3 ,Tcold) — 8 (uj 2 ,T]oot) 

= AJ" — (Thot — Tcold) F.(uJ3, Tcold)- 

In the second line, we employed the thermodynamic iden¬ 
tity 8 — F + TT-L, which follows from the definition of 
F. During the second isothermal stroke, the excess work 
Wex ,3 is performed on the system by the cold reservoir. 
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This allows for the system to exhaust the excess heat 
Qex ,3 at temperature Tcoid- Hence we have 

Wex,3 = -^(W4, Tcold) — -^(^3, ^cold) 

Qex,3 = 7cold('H(aJc, Tcold) — 'H(w3, Tcold)) • 

Finally, during the second isentropic stroke, the cold 
reservoir performs the excess work Wex ,4 on the system. 
No excess heat is exchanged and the temperature of the 
engine increases from Tcom to Thotj 

Wex,4 = AS = £{LOi,Thot) -£(w4,Tcold) 

= at + {That — Tcoid) Thot)- 

As before, Eq. (15) reflects the isentropic condition, 
'H(wi,rhot) = 'H{L04aTao\d)- 

The efficiency of a thermodynamic device is defined as 
the ratio of “output” to “input”. In the present case the 
“output” is the total work performed during each cycle, 
i.e., the total excess work, Wex = W + Qc- There are 
two physically distinct contributions: work in the usual 
sense, W, that can be utilized, e.g., to power external 
devices, and correlation energy, Qc, which cannot serve 
such purposes as it is the thermodynamic cost to main¬ 
tain the non-Gibbsian equilibrium state. Therefore, the 
only thermodynamically consistent definition of the effi¬ 
ciency has to read 



FIG. 2. (color online) Efficiency of the quantum Carnot cycle 
for Eq. (17). The blue, solid line results from W/Qtot, which 
does not properly account for the correlation energy Qc- The 
red, dashed line is the Carnot efficiency (16). Parameters are 
Tcoid = 1, Thot = 1.5, LOi = 0.2, u !2 = 0.6, h = I, and m = 1. 

a. Example 2 This deviation becomes even more 
dramatic in our second example. Consider a quantum 
particle in a quantum piston as in Fig. 1. Such a system 
can be realized, for instance, as a qubit coupled to an 
optical cavity with Hamiltonian {h = 1 ) 

H = -^ az + uJbO^a + g {a^ + a). (18) 


V = 




= 1 - 


Tcoid 

Thot 


= VC, 


(16) 


which is identical to the classical Carnot efficiency. Ear¬ 
lier analyses did not distinguish between correlation and 
excess part of the heat, and the efficiency was simply 
defined as g = WjQtot, see for instance [15]. 

Example 1 To illustrate these concepts and to build 
intuition we now turn to illustrative systems. As a first 
example we consider quantum Brownian motion, i.e., a 
harmonic oscillator coupled to an ensemble of harmonic 
oscillators. In this case the non-Gibbsian equilibrium 
state of a Brownian particle with mass m becomes [27], 


{x\a\y) 


1 

\/ 27 r (a; 2 ) 


exp 


{x + yf _ {x - y)^ \ 
8 (x 2 ) 2hy{p^))^ 

(17) 


where (x^) = {Dpp + m'jDxp )and i^p^) = 
Dpp/^. Here 7 is the damping constant, and Dpp and 
Dxp are diffusion coefficients, which read in a high tem¬ 
perature expansion Dpp = (w^ — 7^)/12 

and Dxp = 13/12 [52]. To implement the cycle we 

assume that the angular frequency uj is controlled exter¬ 
nally. 

In Fig. 2 we plot the resulting efficiency. We observe 
that the Carnot efficiency is, indeed, attained for all val¬ 
ues of 7 if one properly accounts for the correlation heat. 
The blue, solid line is the ratio of work over total heat, 
W/Qtat- Notice that in this case 77 = W/Qtat deviates 
from TIC- 


Here, a, are the annihilation and creation operators 
of bosonic modes with frequency Wf, [53] . The base fre¬ 
quency ujb is the parameter to be changed, which can be 
experimentally realized by varying the laser in the cav¬ 
ity. Pauli matrices Uj, ctz represent a two level atom with 
energy ujq /2 and its coupling to the cavity [4]. Finally, 
the last term ax {a/ + a) = a^, x describes the inter¬ 
action of qubit and piston and can be interpreted as an 
intra-system non-conservative forcing of strength g. 

Finally, the thermal reservoirs are phenomenologically 
modeled by a Lindblad master equation [54], C{p) = 
-i [H, p] + Vq {p) + Vb {p), where 


T, (p) = Iq (iV„, + 1) {^-p(y+ - ^{cr-cr+,p}^ 
^cr+pcr_ - ^{cr+cr_,p}^ , 
with ladder operators for the atom, a±, and 
Tb {p) = lb {Nu,, + 1 ) (^a-pa^ - ^{aa\p}^ 

+ lbNcj, (^a^pa- i{a'fa,p}^ . 


(19) 


( 20 ) 


Here, = l/(exp (/3x) — 1) and iq, ib are fermionic and 
bosonic coupling constants, respectively. Lindblad mas¬ 
ter equations are generally applicable only to describe 
quantum system weakly coupled to the (classical) envi¬ 
ronment. For the present case this assumption is justified 

















5 



FIG. 3. (color online) Efficiency of the quantum Carnot cy¬ 
cle for Eq. (18) with Eqs. (19) and (20). The blue triangles 
are a numerical verification of (16), whereas the green circles 
result from WjQtot, which does not properly account for the 
correlation energy Qc- The inset is a magnification for small 
values of g. Parameters are ojq = 1, 7 ^ = 75 = 0.05, Tcom = 1, 
Thot = 1.5, Lui = 0.2, ijj2 = 0.6. 

as the Hamiltonian (18) describes a generic quantum op¬ 
tomechanical system, for which Lindblad master equa¬ 
tions have been proven to be adequate [54, 55]. More¬ 
over, we do not have to account for dynamical corrections 
as we are only interested in quasistatic, i.e., infinitely 
slow processes. It is worth emphasizing that from micro¬ 
scopic treatment one would expect that the interaction 
between the two subsystems changes the individual dissi- 
pators [56]. However, for a macroscopic Lindblad master 
equation those corrections would force the system to re¬ 
lax into a Gibbs state. For the present purposes, we have 
specifically chosen a phenomenological system which does 
not relax into a Gibbs equilibrium state. 

In this model the two subsystems are coupled to the 
thermal reservoir independently. However, they “feel” 
each other through the direct interaction. Only in the 
limit 5 —> 0 the steady state is a Gibbs state [54]. For fi¬ 


nite interaction qubit and cavity are correlated and they 
share information. The thermodynamic price for main¬ 
taining this correlation during the thermodynamic cycle 
is the correlation energy Qc (8). 

Figure 3 plots the resulting efficiency (16). We ob¬ 
serve again that the classical Garnot efficiency is, indeed, 
attained for all values of g. The green circles are the 
ratios of work over total heat, WjQtot- Notice that in 
this case the Garnot statement appears to be violated as 
r] = WjQtot can be larger or smaller than 77 c as a func¬ 
tion of g. This apparent violation is not physical, but is 
rather rooted in an thermodynamically inconsistent iden¬ 
tification of the excess heat. 

VI. CONCLUDING REMARKS 

The present studied analyzed the thermodynamics of 
non-Gibbsian quantum heat engines - devices that oper¬ 
ate cyclically in non-Gbbsian equilibrium states. We in¬ 
vestigated the thermodynamic processes underlining such 
nanodevices and concluded that it is impossible to har¬ 
ness quantum correlations in quasistatic processes to en¬ 
hance the maximum efficiency of such devices. Instead 
one has to modify the definition of heat, and account 
for the correlation energy necessary to maintain coher¬ 
ence and correlations. In conclusion, we showed that the 
Garnot statement of the second law is universally valid 
also for quantum heat engines. 
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